Equivalence of Classical Skyrmions and Coherent States of Baryons 
II. Baryonic Coherent State Construction on Compact Manifolds 
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In connection with the possibiHty of skyrmion production from small domain disoriented chiral 
condensates formation from heavy ion coIUsions, the direct relation of a classical skyrmion to baryon 
states is examined. It is argued that a skyrmion is a coherent state of baryons. The collective 
coordinate approach of quantization means that the physical baryon states exist not in flat space 
but on a compact manifold. This requires the construction of coherent states in such a curve space. 
Using the techniques associated with the Segal-Bargmann transform also known as the coherent 
state transform used for example in the study of the classical limit of quantum gravity, such states 
can be constructed in the context of the Skyrme model. They are made up directly of baryon states 
on S"^ but with quantum operators on the 5*0(3) manifold. In terms of wavefunctions, they are a 
superposition of the analytic baryon wavefunctions of Adkins, Nappi and Witten. The distribution 
of the baryon states in terms of the relative probabilities of the baryons inside a skyrmion can 
therefore be determined. 

PACS numbers: 12.39.Dc, ll.30.Rd, 25.75.-q, 03.65.-w 
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I. INTRODUCTION 



Recently the possibility of skyrmion formation from 
disoriented chiral condensates (DCC) was raised j|, ||. 
This is in response to the discrepancy of data on the rare 
hyperon yields at the Super Proton Synchrotron (SPS) 
with those from numerical models. The multi-strange 
baryons such as the and f2 also exhibit unusual trend 
when compared to the other hadrons l^i |§7 0, D ■ Com- 
parison with thermal model shows that the deviations 
for Q and S were found to be concentrated in the low 
k± region All these point to an apparently different 
production mechanism of the multi-strange baryons. It 
may very well be that the mechanism is there also for 
the non-strange baryons but they are too abundant for 
the small fractional increase in the yields to be notice- 
able. Skyrmion production is one such baryon produc- 
tion mechanism that can provide an explanation for the 
discrepancy. This channel of generating baryons through 
the Kibble mechanism ||lO| has been in the literature 
for some time ju], yf[l|, |l|l. But unhke for ex- 
ample monopole boundary conditions are important for 
skyrmions and must be properly taken into account . 
An important difference between earlier considerations 
and Ref. ||, ^, |j is that in the latter skyrmion formation 
is connected to that of DCC and to data at the SPS for 
the first time. 

Ever since the first proposal of domain formation dur- 
ing the chiral phase transition could lead to observable 
consequences in the form of fluctuations of_the charged to 
total pion ratios |17, 1^, 19, 
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much exper- 
imental effort has been spent on the search for DCC. No 
sign of them has been found so far |2^, ^ at least 
not via pions. While this may be due to no DCC forma- 
tion, it could also be that only small domains are able 



to form. In the latter case, observation of DCC through 
pion distribution becomes impossible ||2^, This 
is when the possibility of skyrmion formation becomes 
important. Whereas the observation through pions re- 
lies on having large domain size, the opposite is true for 
baryon formation through skyrmion. In the latter case, 
the smaller the domain size, the larger the probability 
of skyrmion formation [p^ . Of course there is a natural 
limit of how small the domain can be. 

In this series of papers, we intend on exploring this 
possibility and try to put this idea on a more concrete 
foundation than Ref. [|], |[ |^ were able to do. One of 
the first requirements is to answer the question of what 
a skyrmion really is besides being a baryon. Evidently 
there are many different types of baryons. If a skyrmion 
is produced, what observable baryon or baryons are ac- 
tually produced? There are various attempts in the past 
to try to project out known baryons from chiral solitons 
for studying properties of the baryons |30[ ^ ^ . 
The chiral solitons themselves are invariably constructed 
out of three valence quarks surrounded by a cloud of co- 
herent pions. In this type of approach, the coherence is in 
the pion cloud and not in the baryon states. Another dis- 
tinct approach is that of Amado et al ^ where aux- 
ihary boson fields are used to reahze the SU (2) x SU (2) 
algebra. From these boson operators, they constructed 
group theoretic coherent states for any given fixed num- 
ber of bosons. This fixed boson number can be identified 
with Nc the number of colors. Thus by going to the 
large Nc limit, these group theoretic coherent states can 
be considered to be a skyrmion. In spite of the success of 
either of these approaches of being able to decompose the 
solitons into known baryons and therefore provide the an- 
swer to the question of what baryons will be produced by 
skyrmion formation, we will use a third approach already 
introduced in Ref. [BTI The chiral soliton approach 
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is based on the linear sigma model and relies on the bo- 
son fields to be linear independent of one another whereas 
the Skyrme model is a non-linear sigma model. Therefore 
the results of the chiral soliton model cannot be applied 
to the Skyrme model. The second approach of Amado 
et al is closer to our method but with an important dif- 
ference. Although both methods are based on coherent 
states, those in Ref. (SSl ^ are of the group theoretic 
type of Perelomov pS|,^[Q| where the states are generated 
by a relevant group action on some fixed state vector. In 
Ref. Adkins et al have shown that the baryon wave- 
functions from the collective coordinate approach are all 
analytic functions on S^. Our goal is to obtain coherent 
states based on these wavefunctions of Ref. |4l). As a 
consequence the suitable coherent states must be those 
on S^. Those of Amado et al although quite flexible are 
not the ones that we are after. We have no desire to in- 
troduce any auxiliar y b oson fields and rather be as close 
to the work of Ref. |Q as possible. 

The coherent states that we will use are not of Perelo- 
mov type but are coherent states on compact manifold 
appropriate to the problem at hand. We have out- 
lined the steps in g^. These coherent states on a 
compact manifold follows from the generalized Segal- 
Bargmann transformation or the coherent state trans- 
form p^ . Their original application were in functional 
analysis and the classical limit of quantum gravity 
|||, m |4|^^||, 1^. But most of the techniques can be 
applied to our present problem. These states constructed 
by using these techniques possess many of the expected 
properties of coherent states. As discussed in ||3^, how- 
ever, straightforward application on will fail and it 
turns out that one must used operators on 5*0(3) but 
states on S^. That is a mixing of ingredient from differ- 
ence spaces is required. These will be elaborated. 

The main result of the paper is that a skyrmion can 
be expressed directly in terms of a coherent superposition 
of baryon states. In the form of superposition of wave- 
functions, these are the analytic functions on derived 
in The probabilities of the distribution of the dif- 

ferent baryons in the superposition can be written down 
analytically in the form of modulus square of wavefunc- 
tion weighed by an exponential factor with the exponent 
given by the ratio of the energy of the particular state 
to a fundamental energy scale. Application of the results 
to actual DCC formation in heavy ion collisions will be 
done in a future work . 

The organization of the paper is as follows. First in 
Sec. II we briefly review the Skyrme model and the re- 
sults derived in [Q. Then we discuss the connection of 
classical solutions to coherent states in Sec. III. The 



main method of constructing coherent states on a com- 
pact manifold will be introduced in Sec. IV and applied 
to SU (2) quantized Skyrme model. As mentioned one 
must overcome an obstacle in applying the method and 
this is discussed in Sec. The solution is given next 
in Sec. VI. Finally in Sec. |VIII , the superposition of 
baryon states will be shown and examples will be given. 



A general expression for the probabilities of the differ- 
ent states will be written down and will be calculated in 
some simple examples. 



II. THE SKYRME MODEL 

The Skyrme model has the Lagrangian density [^9[ |5C 
2 1 



tr[U^d^U,U^d,U]^ (2.1) 



where 



U = exp{iT • 0//^} = {(7 + iT ■ tt)//^ 



(2.2) 



/tt is the pion decay constant and g is the p-ir-n cou- 
pling H. This La grangian density has the classical 
skyrmions as the solutions to the Euler-Lagrange equa- 
tion. 

To quantize the Skyrme model around the maximally 
symmetric classical skyrmion solutions 



U = Us = cxp{ix • f F{r)} 



(2.3) 



time dependence has to be introduced to these static so- 
lutions. One can do this using the well known collective 
coordinate quantization first done by |4ll within the con- 
text of the Skyrme model. There are two ways to do this. 
The well known one is by using SU{2) collective coordi- 
nates 



A = aa + ia ■ T . 



(2.4) 



Here A is an element of SU{2) where the components are 
functions of time Ub — ab{t) and subject to the unitary 
constraint 

Co + - 1 = abflfe -1 = 0. (2.5) 

Because of the form of Eq. ( ^.l| ) given a solution Us, 

U's^AUsA^ (2.6) 

is also a solution. In terms of A the Lagrangian takes the 
form 



L = 2\abab - M 



(2.7) 



M is the basic mass scale of the skyrmions and A has the 
dimension of inverse mass. The explicit integral expres- 
sions of these constants in terms of F(r) can be found for 
example in Ref. [^l| [5^ and is also given in our previous 
paper p8[ . The conjugate momentum to the coordinate 
Of, is naturally given by Hb — 4Adb. The Hamiltonian is 
therefore 



H 



1 

8A 



TTbTTb + M 



(2.8) 
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Proper quantization using the Dirac bracket 1 53 , 5^ gives 



[afc,ac] = (2.9a) 
[ab,7rc] = i{Sbc-abac) (2.9b) 
[TTb,TTc] = i{dcTTb - dbT^c) (2.9c) 



Another approach is to map from SU{2) to 5*0(3) and 
quantize using the S0{3) coUective coordinates. This ap- 
proach is known but as far as we are aware, it has never 
been done formally or shown explicitly. In the literatures, 
one can at most see a piece of information here and an- 
other one there. In we have shown how to do this 
rigorously using the Dirac quantization procedure when 
dealing with a constrained system. In this case, one uses 
the map 



At, At = TjRji{A) 



(2.10) 



to map an element A of SU{2) to an element R of SO{3). 
The constraints in this case are the requirement of unit 
determinant 



det i? = 1 



(2.11) 



the identity is also an element of the group and the in- 
verse is given by the transposed matrix 



RR- 



R-^R = 1 and R^^ = R^ 



(2.12) 



After some algebraic manipulation, one can arrive at the 
Lagrangian 



L — jX RijRij 



M. (2.13) 
With the conjugate momenta to i?y evidently given by 

(2.14) 

the Hamiltonian follows easily 

M . (2.15) 



ij - i ny iiy 



The quantization using S'0(3) collective coordinates yield 
some unfamiliar but interesting commutators 



[Rij,Rki\ — 



(2.16a) 



[Rij.tlkl] = \i{5ik5jl + eikm^jhiRmn — RijRkl) 



[iiy,nfci] — ji5ik{Rml^mj — Rmj^ml) 
~^'^^^jl(^Rkm^iTn Rim^ki^ 



(2.16b) 



(2.16c) 



Details of these derivations have been laid down in . 

It has been shown in that in either case, the clas- 
sical Hamiltonian can be expressed universally as 



H 



1 

2A 



J' 



M 



1 

2X' 



M 



(2.17) 



independently of whether SU{2) or 5*0(3) collective co- 
ordinates were used for the quantization. Quantization 
in both cases result in replacing J and / in Eq. (2.17) 
by their operator counterparts. The differential form 
of the operators are however quite different |Q. The 
same is true for the resulting eigenstates. In the case 
of using SU{2) collective coordinates, quantization pro- 
duce states with both integral and half-integral spin and 
isospin states. Their wavefunctions belong to L'^{S^) and 
are polynomials in {ab + iac). The energy of these states 
depend on the degree of the polynomial I via 



El = + 2) + M . 



(2.18) 



Only the half-integral spin and isospin states can be inter- 
preted as physical baryon states. The lowest lying states 
are shown in Appendix Quantization with 50(3) col- 
lective coordinates, on the other hand, produces only in- 
tegral spin and isospin states. None of these are physi- 
cal particles [Q. The wavefunctions are polynomials of 
products of combinations of Again the energy of 

the states depends on the polynomial degree I but with 
a slightly different dependence 



El = l^l{l + l) + M . 



(2.19) 



III. 



CLASSICAL SOLUTIONS AND COHERENT 
STATES 



Although the static skyrm ion solution w as us ed to ob- 

and Eq. ( 2.15 ), its only 



tain the Hamiltonian in Eq. (2 



presence is in the two constants A and M . The clas- 
sical solution determines the mass spectrum Ei of the 
quantum states but does not provide a connection of 
the classical solutions to the latter. To find a connec- 
tion, recall in the case of the simple harmonic oscillator 
(SHO) that Schrodinger showed a long time ago that the 
quantum analogs to the classical solutions can be found 
Ipql - These are the so-called coherent states which are 
superposition of the eigenstates of the quantum system. 
Therefore our present problem should be to find coher- 
ent states of baryons which are analogs of the classical 
skyrmions. If this can be done then whenever a skyrmion 
is formed, the quantum analog will provide information 
about the distribution of baryons that can be produced 
via DCC. The question is how to do this. 

There exist in the literature many papers that touched 
upon coherent states, but most of which are for applica- 
tions far removed from our problem (see for example Ref. 

This just shows how coherent states are relevant to 
many different areas of physics. Only those that gener- 
alize the concept of coherent states outside the context 
of the SHO are potentially relevant. The closest works 
in the literature are the chiral soliton model and the in- 
teracting boson model treatments already mentioned in 
the introduction. There are at least three principal direc- 
tions for generalization: one that is based on the action 
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of a group, another that is based on minimizing the un- 
certainty relation and a third is based on eigenvectors of 
the annihilation operators. 

The main proponent of the group based method is 
probably Perelomov |39, ^ but there are others (for 
example see the references in Ref. This general- 

ization is based on that aspect of the original coherent 
states where they can be considered to be generated by 
the action of the Heisenberg-Weyl group on the vacuum 
vector. That is 

IV') = D{a)\0) = exp(ait „ a*A)\0) . (3.1) 

where A and A' are the annihilation and creation opera- 
tors respectively. Apart from an additional trivial phase 
factor, D{a) is a representation of the Heisenberg-Weyl 
group. The generalization in this group based method is 
then 



(3.2) 



where g is an element of a group G with a representation 
T{g)- \ipo) is some fixed vector on which the coherent 
states are constructed. This generalization is not par- 
ticularly suited to our problem since a \ipo) has to be 
picked which can be quite arbitrary. Also in the case of 
G = SU{2) X SU (2), which is relevant to our problem be- 
cause of spin and isospin, the action of the group on the 
states will only be restricted to a multiplet of a particular 
total (iso)spin j if l-^o) is in the j multiplet. If iV'o) on 
the other hand is already a superposition of a number of 
state vectors from different multiplets, the arbitrariness 
of IV'o) will be further increased. The fact of having to 
choose a fixed iV'o) at the start already affects to a large 
degree what the coherent state will be. In our problem, 
this is tantamount to selecting some baryon states and 
force the skyrmion identity on them. This is undesirable 
and there is little justification for that. 

The second approach based on minimizing the uncer- 
tainty relation to AxAp — h/2. The resulting coherent 
states can then be labeled by the ratio Ax/Ap since their 
product is fixed at the minimum | |59| . However this ap- 
proach is not terribly appealing in the sense that the an- 
nihilation and creation operators depend on energy lev- 
els. Also why should the alternative, say AxAp = h or 
AxAp — 3h/2, would be bad choice for coherent states 
psf . This approach was proposed to deal with different 
type of possible potentials. It is not clear how one can 
implement this on compact spaces. In the last approach 
coherent states are simultaneous eigenstates of the an- 
nihilation operators. As criticized in Ref. there is 
little value from the physical point of view why such a 
mathematical criterion alone should automatically give 
the resulting coherent states physical meanings. 

The above discussion made clear that relying on only 
one aspect of the original coherent states and generalizing 
based on that alone is not a physically sound procedure. 
Much stricter criteria must be imposed which require, at 
the least, more than one of the above aspects of coherent 



states to hold. It should be reasonable to find coherent 
states with properties close to or are the generalization 
of those of the simple harmonic oscillator. This need 
not be true in general because there are many different 
quantum systems that bear little resemblance to that of 
the SHO. Only some features of the original should carry 
through to these other systems. I n ou r case the mani- 
fold is S^. The Hamiltonian in Eq. (2_^) has no potential 



so classical motions are geodesies or great circles on the 
four-dimensional sphere. As is well known that the one- 
dimensional projection of this motion onto any of the Ub 
axes is SHO-like. The only difference here is the motion 
in the different dimensions are not completely indepen- 
dent. One can see the similarity in the Lagrangian if 
the constraint Eq. (2.5) is included using the Lagrange 
multiplier method 



L = 2Adfcdb — p{a(,ab — 1) — M , 



(3.3) 



where p is the Lagrange multiplier. Apart from the con- 
stant terms, this resembles the four-dimensional SHO La- 
grangian. It should therefore be reasonable for the orig- 
inal coherent states of the SHO to be a guide for the 
construction of coherent state in later sections. 



IV. CONSTRUCTING COHERENT STATES ON 
A COMPACT MANIFOLD 

From the discussions in the previous sections, it is ev- 
ident that coherent states satisfying the right physical 
criteria are required for the problem at hand. Addition- 
ally unlike most other situations the coherent states must 
be on a compact manifold of S^. Fortunately there is a 
method to do this. For example it has been done on a 
circle and a two-dimensional sphere by Kowalski et al 
p3| , p4[ , and by Hall et al using another approach on 
general S*" Kowalski et al based their method on 

the relation between the polar coordinates and the an- 
gular momenta, and the analogous relation between the 
conjugate pair {x,p) in the SHO. In the standard case 



X + ip 



so we have 



and 



a\a} = a a) 



p/3(S+jp) _ p/3q 

The analog of this on a circle is 



(4.1) 



(4.2) 



(4.3) 



(4.4) 



Kowalski et al method exploited this analogy. Hall et 
al, on the other hand, used techniques developed in the 
studies of analytic functions on compact manifolds and 
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applied to canonical quantum gravity, namely the gen- 
eralized coherent state transform [|42[ together with the 
complexifier method of Thiemann [ [46| . We will use this 
second method below since it is much easier to generalize 
to any dimensions. In spite of the different approaches 
and from completely different points of view, the results 
of Ref. turned out to be two special cases of those 

in Ref. [p5| . This strengthens the method and provides 
some significance to the so constructed states. 



A. The Method 

According to these works, one starts by constructing 
the annihilation operators. U X = (Xi,X2, ■ ■ • ) are co- 
ordinates on a compact manifold then the annihilation 
operator A — (^1,^2,...) and the creation operators 
are given by 



C 



Ai = WX^W-^ = e-^X^e 
ij = W-'^X,W = e'^X.e-'^' 



(4.5) 
(4.6) 



where W — e~ is there for phase space Wick rotation 
(in quantum gravity) with C the complexifier operator 
being the generator of this transformation. Given the 
annihilation operators, simultaneous eigenstates of these 
operators can be determined and the coherent states will 
follow from them. 

Eq. ( |4.5| ) has its origin from the classical version of the 
phase space Wick rotation. Classically the complexifier 
function C (to be defined below) generates the trans- 
formation from the real phase space pair {X, P) to the 
complex {X"'' , P"") pair via 



n=0 



(4.7) 



and similarly for P^. Here 

{X, C}n+1 = {{X, C}n, C} = {... {{X, C}, C},...,C} 

v ' 

n 

(4.8) 

with {X, C}o = X. The latter denotes applying the Pois- 
son bracket n+l times. Quantization promotes functions 
to operators and Poisson brackets to commutators 

0^0, {01,02} ^ -[61,62]. (4.9) 

Therefore 



X^ = J2±[X,CU = e-^Xe' 



(4.10) 



n=0 



with [X,C]o = X and [X,C]„+i = [[X,C]„,C] repre- 
senting the repeated application of the commutator n+l 
times. The complexified X is the classical annihilation 



function X'^ = A. Eq. (4.5) is therefore the quantized 
version of the transformation. 

From Ref. the complexifier function C itself is 

given by the kinetic energy term in the Hamiltonian di- 
vided by an energy scale so that C is dimensionless. This 
scale defines a fundamental energy scale to the problem 
under consideration ]4^ . Introducing u as this scale and 
taking the kinetic energy from a general Hamiltonian de- 
composed into kinetic and potential terms 



H = T + V , 

the complexifier function is 

T ^ f 
C = >C=- 



(4.11) 



(4.12) 



which becomes the operator C upon quantization. 

The annihilation operators are now completely defined. 
The coherent states are found by first finding simultane- 
ous eigenstates of these operators. On a manifold, one 
can always define the position states \x) by 



X,\x) = 



(4.13) 



where x = {xi,X2, ■ ■ ■) represents a point on the mani- 
fold. From the form of Eq. (4.5) simultaneous eigenstates 
of Ai can be constructed from |a;). This is done by acting 
on \x) with the Wick rotator W 



W\x) = e-^\x) (4.14) 



so that 

A^ 



= A^W\x) = WX^\x) = x^W\x) = x^\^) . (4.15) 

are now simultaneous eigenstates of Ai. They are, 
however, not yet the coherent states that we are looking 
for. Remembering from the SHO, the eigenvalues of the 
annihilation operators are in general complex but the Xi 's 
are real. The complex nature of the eigenvalues contain 
information both on positions and momenta. On the 
other hand eigenvalues Xi being real can carry only half 
the information. To obtain the coherent states from 
one must perform an analytic continuation from real x 
to complex x^ . Once this last step is taken, the coherent 
states are 



W\x^ 



with 



(4.16) 



(4.17) 



The eigenvalues in the set x^ = {x^, x^, . . .) serves also 
as the label of the coherent states. It follows that there 
is one coherent state for every point in phase space, the 
same as in the SHO. 

Superficially we are following the third direction for 
generalization of the SHO coherent states as discussed in 
Sec. [II. In reality much more stringent conditions come 
automatically with this method. This has been addressed 
in a couple of papers by Thiemann et al ^ . 
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B. An Illustration with the Simple Harmonic 
Oscillator 



The last subsection described the method of how to 
arrive at the coherent states for a given manifold. Let us 
illustrate it with a familiar example of the n-dimensional 
SHO 0. In t his c ase the complexifier is C = 

/2mu). From Eq. (4.10) with the usual commutation 
relations 



i6, 



the i component of the annihilation operator is 

i 



A,. ^ XI 



X,, 



muj 



(4.18) 



(4.19) 



Up to a factor of muj / 2 these are the usual annihilation 
operators. According to the previous section, a coherent 
state labeled by the complex vector r"' = (rf , , . . . , r^) 
is 



(4.20) 



In this form, it does not resemble a coherent states of 
the SHO. To recover the familiar form, one goes to the 
position representation 



Inserting now a complete set of momentum states 



= / d> {x\e- 



(27r) 



n/2 



(4.22) 



Writing = r + ip, the wavefunction becomes 



V 27r / 

(4.23) 

After normalization this is the usual coherent state wave- 
function 



(a;|'0, 



n/4 



^imujp-{x — r) ^—muj{x — r) /2 



It has the n-dimensional vector of eigenvalues 



(4.24) 



v/mw/2rf = ^JnlUJ/2{r,+ipi) (4.25) 

for the correctly normalized annihilation operators. Ob- 
serve now that the introduced fundamental energy scale 
m is the fundamental frequency (up to K) of the oscillator. 

In the case of the SHO, the properties of the coherent 
states are: 

• there is one coherent state per phase space point, 



• they saturate the Heisenberg uncertainty bound, 

• the wavefunction peaks at the location given by the 
label in the classical limit, and 

• Ehrenfest theorem is satisfied. 

The first property was seen in the previous section. 
The remainders have been painstakingly shown to hold 
for general coherent states on compact manifolds con- 
structed from this method in the classical limit in Ref. 
pO| , |6H . The actual proofs are lengthy and involve rather 
heavy mathematics. 



V. 



ARE COHERENT STATES ON 
SKYRMIONS? 



Returning to the problem at hand, from the Hamilto- 
nian in Eq. (^.8[) the complexifier is 



(5.1) 



In terms of the position operators on which are now 
the fih, the annihilation operators are 



Ob 62 



(5.2) 



The form of Af, does not resemble Eq. ( 4. 19] ) at all but it 
can be brought into a form much closer in resembl ance . 
This is shown in Appendix O. A s seen there Eq. _([5.2| ) 
is really the analog of Eq. (^.19| ). In Appendix ^ we 
address the uniqueness issue of Ai,. The corresponding 
eigenstates labeled by complex variables are 



(5.3) 



According to the method, this is the general form of the 
coherent states on for our Hamiltonian. Mathemati- 
cally the problem is solved. Unfortunately on the physi- 
cal level, we have not reached our goal at all. 



A. The Problem with the Coherent States on S 

As discussed in Ref. |3^ we would like to expand the 
skyrmion in terms of baryon and excited baryon states. 
Eq. ( ^.3| ) is not yet in this form. As mentioned in Sec. II 
the SU (2) quantized Skyrme Hamiltonian admits both 
integral and half-integral (iso)spin states and not a ll o f 
them are physical. One has to first verify that Eq. (5.3) 
carries only half-integral (iso)spin states before it can be 
identified with a classical skyrmion. In Ref. it was 
shown that applying the method on failed to yield a 
quantum analog to the skyrmion. It fails precisely due 
to the presence of unphysical integral (iso)spin states. 

The complete Hilbert space on which the Skyrme 
Hamiltonian acts is spanned by |j, m, n) where j, m,n € 



7 



a, 



a-. 




the states |j, m,n). In the Appendix g some algebra of 
Ji and li with a;, have been worked out. It is convenient 
to use the combination 



CLi ± 10,2 



(5.6) 



The algebra shows that they raise or lower the third com- 
ponent of spin and isospin of \j, to, n) by 1/2. For exam- 
ple 



FIG. 1: This figure summarizes the and . Each of 
which is a member of simultaneous spin and isospin j — 1/2 
multiplet. J- will bring at down to and /+ moves aj up 
to a7 etc. 



haflj, 771,71) 



(m± i) af\j,7n,7i) 



(5.7) 



T 2) «t b';'71,7l) (5.8) 



Z.Z- 



but the d\ j, TO, n )'s a re not eigenstates of and 1^ in 
general (see Eq. jB9[ )). The algebra can succinctly sum- 
marized by Fig. This shows that each of the and 



1/2. Inserting a complete set of states in Eq. (|5J) g± belongs to a simuhaneous spin-1/2 and isospin-1/2 



IV', a^) = H 2L-f'|j,TO,7l)(j,771,77|a'^) 

j,m,TieZ,Z+l/2 

E e-^^-(^-+i)^;„„(aC)b-,TO,7i) 



j,m,neZ,Z+l/2 



(5.4) 



Here (j, TO,7i|a^) = ''/'jmn(o^) the complex conjugate 
of the wavefunction of the state with quantum numbers 
(j, 771, 77) evaluated at the point on the complex four- 
dimensional sphere (the complex conjugation is per- 
formed first for real a before a = a*' is set). A complete 
set for j = 0, 1/2, 1, 3/2 were given in Ref. |Q and only 
the j = 1/2, 3/2 ones are repeated in Appendix As 
discussed in [ p8[ and in Sec. II all wavefunctions are 
polynomials made up of {ai, +iac) on 5*^. One can there- 
fore analyticall y co ntinue them to the complex S^^. The 
expansion Eq. ( ^.4[ ) shows manifestly the presence of the 
unphysical states. We can conclude that \tp,a''^-^) is defi- 
nitely not the classical analog of a skyrmion. 



B. Projecting Out the Physical States? 

Can the situation be salvaged somehow? How a bou t 
simply dropping all the unphysical states from Eq. (5.4) 
or imagine that we have a projection operator V at our 
disposal, which when applied to \ip, a'^), projects out only 
the physical states 



j,m,neZ+l/2 



\j,7n,77) {j,m,n\a^ 



(5.5) 



In Ref. |37| an argument in terms of the wavefunction 
representation as to why this would not work was pre- 
sented. Here we shall instead rely on the algebra of S^. 
The form of Eq. (4.5) requires us to know how ai, acts on 



doublets. 

Their actions on the state vectors in general are 



a^|j, 777 , 77) 



at\j,7n,7i) 



CtU,77l,7l,+)\j + i,TO± i,77± i) 
+ Cf{j,7n,7l, - i,TO± 1,77± 1) 

(5.9) 

Ct{j,77l,7l,+)\j + i,TO± i,77=F I) 
+ Ct{j,77l,7l,~)\j ~ i,TO± i,77=F I) 

(5.10) 

where Cf^ and are numerical coefficients and func- 
tions of {j,7n,n). These equations can be formally de- 
duced from the algebra of the SU (2) collectively quan- 
tized operators. These are given in Appendix |b| How- 
ever the simplest way to see this is to use the observation 
that each af and af is a member of a doublet. Therefore 
they have j = 1/2 and to = 77 = ±1/2. The eigenstate on 
which the operators act also has its own value of j, to, 77. 
The operators and the state thus form a direct product 
which can be decomposed as a direct sum of state with 
new total spin and isospin eigenvalues j -t- 1/2 and j — 1/2 



J (g> 1/2 = [j + 1/2) © (j - 1/2) 



(5.11) 



This results naturally in Eqs. (|5.9| ) and (|5.10|) . In Ap- 
pendix ^ more details are given. Also given there are 
the necessary steps to solve for the coefficients and 

Acting on Eq. (q^) with ai, is to turn every half- 
integral spin state in th e ex pansion into integral spin 
state. Therefore Eq. (5.5) is not an eigenstate of 



Eq. (5.2). Keeping only physical states in the expansion 
does not solve the problem. 



VI. COHERENT STATES WITH 50(3) 
OPERATORS ON THE S'^ MANIFOLD 

We have seen that the problem of constructing coher- 
ent state using the SU (2) collective quantization and that 
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the 5*0(3) quantization giving only unphysical states. 
Thus neither approach permits the identification of the 
coherent states as the quantum analog of the classical 
skyrmions. Nevertheless there is a way out. As discussed 
in Ref. |^ the reason that the SU{2) theory fails was be- 
cause of the SU{2) operators mapped fermions to bosons 
and vice versa. Therefore the SU (2) annihilation opera- 
tors require the full set of bosonic and fermionic states. 
However if we act on the Hilbert space of the SU{2) quan- 
tized theory exclusively with combinations of a^aj, this 
would map fermions to fermions and bosons to bosons. 
The fermion part and the boson part of the Hilbert space 
therefore decouple and the latter can be eliminated. This 
is achieved by "mixing" the SU{2) and 5*0(3) theory. 
One discards the 50(3) Hilbert space but keeps the oper- 
ators and at the same time keeps the Hilbert space (only 
the fermion half of it) of the SU (2) theory but discards 
the operators. 



According to the description in Sec. IV A after intro- 
ducing an energy scale uj, from Eq. ( ^.17 ) the complexifier 
for the 5*0(3) system is 



C* 



i-n,n. 



1 

2Xuj 



J' 



(6.1) 



The annihilation operators based on the 50(3) coordi- 
nate operators Rij are therefore 



Aij — e 



(6.2) 



Just like in the SU (2) quantization, the Aij can be ex- 
pressed in another form by expanding the exponential 
and then regrouping the terms. This will bring them 
to a more traditional form. Using the 50(3) algebra in 
Appendix ^, this other form is derived in Appendix |^. 

The coherent states are constructed from a complex 
labeled position state |a^) on 5*^ as before. This should 
be an eigenstate of the position operators Rij . How do 
the Rij act on \a^)7 Recall that the Rij although are 
elements in the 3x3 rotation matrix R, they are not 
completely dissociated from the at of the SU{2) element 
A. In fact as shown in , Rij = Rij (A) = Rij (a) so 



R,,\a^) ^ R,,{a^)\a^) 



(6.3) 



The position states |a*') are indeed eigenstates of Rij. 
Therefore the method together with the above discussion 



give us 



as the coherent states. So with Eq. 
tion operators we have 



R^ 



')|*,a' 



(6.4) 

as the annihila- 
) ■ (6.5) 



To obtain these states in terms of the eigenstates \ j, m, n), 
one insert a "comple te" set of states in the fermionic 
space to recover Eq. (5^). So finally the coherent states 
consist only of baryon and excited baryon states that we 
have been looking for and identifiable with the classical 
skyrmions. 




FIG. 2: This figure summarizes at least two important facts. 
The first being the action of Ji and li on the nine combi- 
nations of Rij. They form three overlapping spin and three 
isospin triplets. For example J"*" brings R^" horizontally to 
R+° and /" brings R ^ vertically down to R " etc (see alge- 



bra in Appendix Eq. ( F4 ) . Each member is simultaneously in 
one spin and one isospin multiplet. The second is how these 
combinations act on the states (see text). For example i?"*"^ 
raises both the third component of spin m and isospin n by 
one, i?"" lowers m by one but leaves n alone. 



VII. A MIXED QUANTUM SYSTEM OF 5*0(3) 
OPERATORS AND SU{2) HILBERT SPACE 

To permit ourselves to verify explicitly the eigenvalue 
equations and to familiarize with this mixed operator- 
state system, we will derive the action of the 5*0(3) op- 
erators on the SU{2) eigenstates. This can be rigorously 
deduced from the algebra of the 5*0(3) operators. They 
are given in the Appendix ^. From the algebra, again 
one can find useful combinations of the components Rij . 
These have been given in |3^ and repeated in operator 
form in Appendix ^. Each RP'' with p, g = ±, is a sum 
of components Rij and belongs simultaneously to one 
spin and one isospin triplet (see Fig. |^). Using similar 
argument as before, the action of R^'^ on an eigenstate 
|j, TO, n) results in a direct product between the operator 
and the state which can be decomposed as 



J ® 1 = (i + 1) ® j e (j - 1) 

Therefore we must have 



(7.1) 



R'"'\j,m,n) ^ C'"'{j,m,n,+)\j + l,m+pl,n + ql) 
+CP'^ij,m,n,0)\j,m+pl,n + ql) 
+CP'^{j,m,n,~)\j -l,m+pl,n + ql) . 

(7.2a) 

The pi, ql objects need some explanations. The p, q are 
the superscript of i?^* and are signs equal to -f , — or 0. 
pi means the sign p multiplying with 1 . If p is equal to 
-f-, — or 0, then pi = +1, —1 or respectively. The same 
applies to ql. The coefficients are solved and listed 
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in Appendix |^ and | respectively. The actions of 
on |j, m, n) can be similarly expressed but with different 
coefficients. 

Expectation values are an integral part of any quantum 
system and it is most convenient of work in terms of wave- 
functions. The position wavefunction of the coherent 
state labeled by a"" on at the point b ^ (60, 61, 62, ^3) 
on would be (up to normalization) 

*„c(6)-(fe|vI/,aC) . (7.3) 

We use Rij for position operators so the expectation 
value of "position" would be 

(*,aC|i?y |*,aC) = J d%dil~bl) %(6)|vl/,c(6)p . 

(7.4) 

This might cause some confusion because the wavefunc- 
tions are on while the expectation values are of the 
S0(3) operators. In general expectation values of an op- 
erator O from the 50(3) quantized theory in a coherent 
state is 

(vI/,aC|o|v|/,aC) ^ J d%S{l~bl)M>,.ibrd^Ab) ■ 

(7.5) 



VIII. CLASSICAL SKYRMION AS A 
COHERENT SUPERPOSITION OF BARYON 
AND HIGHER RESONANCE STATES 

Given that the coherent states are of the form 

\^,a^) ^Ve~^J'\a^) , (8.1) 

inserting a complete set of fermion states and they can 
be written in terms only of baryon states as 

|*,a^) = J2 e-^^^^+^^\j,m,n){j,m,n\a'^) 

j,m,neZ+l/2 

= E e-^^-(^-+i)^;_(aC)b-m,n). 

j,m,neZ+l/2 

(8.2) 

V'jmn(a) again are the baryon wavefunctions on and 
i^jmni'^'^) the complex conjugates of these wavefunc- 
tions evaluated at the complex point a"". 

Let us see some explicit examples of what the states 
look like. Using the wavefunctions for p,n (j — 1/2) 
and A (j = 3/2), the baryonic coherent state with label 
a^ = f]= (0,1,0,0) is 



= ^(|p,^)-|n,-|))+^^^^(|A++,|)-|A+,i)-K|A",-i)-|A-,-|))+... . (8.3) 

The ellipses denote half-integral spin higher states. In this case the nucleons have the same probability and this is 
similarly true among the A in the superposition. The spin and isospin are correlated by m = n in this special case. 

With a different label = (70, 0, 0, 73) where 70 and 73 are in general complex numbers satisfying 7^^-72 = 1, 

the state is 

■ 3 

l*,7> = — - — ((7o+i73)b,-5> -(70-173)1"-, I) j 

- '^"^ ((70 + «73)'|A++, -I) - (70 + Z73)|A+, -1) (70 - Z73)|A°, \) - (70 - Z73)'|A-, §)) 

+ ... . (8.4) 



Recalling that 7^ are complex so this time both the nu- 
cleons and the deltas do not have the same probabilities. 
This is true in general. However the special choice of the 
label 7 ensures that this is a different correlation m — —n 
in the superposition. 

It should be clear that the expectation values of these 
states depend in general on the label. For example 

(vI',/3|J,|*,/?) = (vl/,/3|/,|vl/,;3)=0 (8.5) 

but 

(4',7|J3|^',7> = (^',7l^3|*,7> ^ . (8.6) 



In the chiral bag soliton type model |3^, ^ |3^, ^ij, 
the construction often involves the restriction of "grand 
spin" to 

(J,+/,)|V)=0 (8.7) 

for minimizing the energy of the hedgehog state [^ . 
From the last expression this can, if need be, be imposed 
by setting 71 = 72 = 73 = and 70 = 1. 

To determine the probability of baryons being pro- 
duced from skyrmion in DCC formation, one can calcu- 
late the relative probability for a given baryonic coherent 
state from the general expression for the states above. 
The relative weights between nucleons and Deltas in a 
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state with a given label is therefore given by 



J2m,n V'l*/2,m,n(a^*)V'l/2,r„,«(a'^ 



The relative probabilities of other higher states can be 
worked out in a similar manner. In the first example 
with ~ 13, the relative probability between deltas and 



nucleons is particularly simple 



= 4 e" 



(8.9) 



In the second case where a = 7, let us take 
Ini7o =Re73 — 0, Re7o — /i and Im73 — v — jj? — 1 
then 



J 



^a(7) 



3_ / (/i + + vf + (/I + vf + (Ai - vf 



I 



(8.10) 



This shows that in general the expression can be quite 
complicated but mostly one would be interested in the 
numerical values and these can be worked out from the 
formula provided A and w are fixed. The former can be 
estimated from the mass of the nucleon, of the delta and 
the pion decay constant. The latter is not known but a 
rough guess would be the energy scale of the chiral phase 
transition or the value of the temperature T^- 



IX. SUMMARY AND OUTLOOK 

We have argued that the classical skyrmion solutions 
can be identified with their quantum analogs, namely 
coherent states of baryons. Due to the non-linear na- 
ture of the Skyrme model, the space of the states is the 
curved, compact space of . Using the method of Ref. 
H, H, H, H, specially suitable for compact 
manifolds, these special superposition states of baryons 
have been successfully constructed. The states and wave- 
functions are exactly those derived by Adkin et al |^ 
without modification including the unphysical states of 
integral spin and isospin ]62| . Since a skyrmion must be 
made up only of baryons, such states must completely 
be removed from the superposition that ultimately will 
be equated to the skyrmion. In order to overcome this 
problem, we find that it is necessary to bring in the opera- 
tors from the lesser known S'0(3) collective quantization. 
Only with these integral spin and isospin operators can 
the fermionic and the bosonic part of the Hilbert space 
be decoupled. The unphysical states can therefore be 
discarded. The distribution of the baryon states of given 
quantum numbers in the superposition come in the form 
of the moduli square of the baryon wavefunction of the 
corresponding states weighed by an exponential factor. 
This factor depends on the ratio of the energy of the in- 
dividual state and an energy scale fundamental to the 
problem. This scale have not been determined but we 
expect it to be of the value of the chiral phase transition 
temperature. With this successful completion of the de- 
composition of the skyrmion into know baryon states, we 



should be ready to apply this to study skyrmion forma- 
tion from DCC in heavy ion collisions. This will be done 
in the near future Eal . 
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APPENDIX A: NUCLEON AND DELTA 
WAVEFUNCTIONS 

For nucleons, 



(ab,+l) 


1 , 

= +-(ai 

TT 


+ ia-}) 




= (flo 

TT 


- ias) 


(a|n,+i) 
(a|n,-i) 


= +-(ao 

TT 
1. 

= (fll 

TT 


+ ias) 
- ma) 



For Deltas, 



V2, 



(a|A++,+|) = (ai+zas) 

TT 



(Al) 
(A2) 
(A3) 
(A4) 

(A5) 



HA+,+i) 



(a|A-,+i) 



i — {ai + ia2Y {ao + ias) (A6) 

TT 

(oi + ia2)(ao + iaa)^ (A7) 

TT 



. ,3 

-I (ao + las) 

TT 



(A8) 
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/fi 

(a|A++, +1) = -i (ai + ia2f{ao - zas) (A9) 

TT 



/2 

a| A+, +1) =-^{a,+ m2)(l - 3(0^ + ^2)) 



(AlO) 



.V2 



(a|A^ +1) = i — iao + m3)(l - + ai)) 

TT 

(All) 

(a|A-, +1) = (ai - m2)(ao + ias)^ (A12) 

TT 



(a|A++,-i) = 
HA+,-i) = 

(a|A",-i) = 
(«|A-,-i) = 

(a|A++,-|; 

(«|A+,-i; 

(«|A°,-i: 

HA-,-f; 



V6 



(fli + m2)(ao - ifls)^ (A13) 



-i (ao - ia3)(l - 3(0^ + al)) 

TT 



V2 



(A14) 

(ai-m2)(l-3(ao + al)) (A15) 



i — (ai — ia2)^(ao + ias) (A16) 

TT 



■V2 . 3 

I — (ao - tas) 

TT 



(A17) 



= (ao - W3)^(ai - m2) (A18) 

TT 

= -z — (01-102)^(00-103) (A19) 

TT 



= (oi - ia2) 

TT 



(A20) 



APPENDIX B: THE ALGEBRA ON S'^ 



In the t ext the action of Ji and li on 05 are required. 
From Eq. (2.9) the commutator between the rotation gen- 
erators (angular momenta) 



Lbc — abTTc — ClcTTb 

and Ob can be deduced 

[Lbc,a.d] = iiSbdac - Sdab) 



(Bl) 
(B2) 



It foUows that the spin and isospin operators and ab sat- 
isfy 



[Ji,ao] = ^icLi , [h,aQ] = 
[Ji,ai\ = -^iao , [li, flj] = 



-lao 



(B3) 
(B4) 



[Jiytt-j] = \ieijkak , [ii,aj] = -ieijuak ■ (B5) 



These can be written more compactly as 

[Ji,ab] = ^i{Sboai - Sib&o + Sbje^jkiik) (B6a) 
[Ii,ab] = ^i{5ibao - SboUt + Sbje^jkcik) ■ (B6b) 

In this form, they are not particularly useful. Let us 
write instead 



fli ± ?a2 



ao ± 10.3 



(B7) 



These commute with each other because the ab's do. 
Their action with respect to spin and isospin are dictated 
and made clear by the following algebra 



[J±, 




= 0, 


[4 




= 


(B8a) 


[J±, 


«?] 


= -ia? , 


[4 


&J] 


= iaf 


(B8b) 


[J±, 


af] 


= 0, 


[4 


af] 


= 


(BBc) 


[J±, 




= iaf , 


[4, 


at] 




(B8d) 


[J3, 






[4 




2 " 


(B8e) 


[J3, 




= H^t 


[4 




1~± 


(B8f) 



The last two lines show that af and af act like raising 
and lowering operators of one-half instead of the usual 
one with respect to the third component of spin and 
isospin. But states acted on by af and af are not eigen- 
states of J2 and P in general. The commutators are 







= -tafJ±±af{J3±l) 


(B9a) 




af] 


= tafJ±±afiJ3±^) 


(B9b) 


[I' 


af] 


= tafi±±af{h±^) 


(B9c) 


[I' 


af] 


3 

= -iafhTafiisT^) ■ 


(B9d) 



Similarly one can work out the commutators with the 
conjugate momenta itb to find out how they act on the 
states. Beginning with the four-dimensional rotation 



[Lbc, TTd] = iiSbdT^c — ScdT^b) ■ 



(BIO) 



This is exactly Eq. (B2) but with db replaced by ifb- 
Therefore itb must have the similar commutators with 
J±, J3, I± and /3 as ab. 



APPENDIX C: ACTION OF at ON THE SPIN 
AND ISOSPIN STATES 

In constructing coherent states on 5*^, it is necessary to 
find out how the operators ab act on the spin and isospin 
states \j,m,n). These are largely governed by the alge- 
bras given in Appendix ^ and the constraint Eq. ( p.5| ). 
Naturally the operator form of the constraint 



abOb 



aja, +dlat 



1 



(CI) 
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is the Casimir operator. The action of Ob on the states 
must take that into account. It must also constrain the 
form that this action will take. From the algebra, one 
can see that they either increase or decrease the m and n 
by one-half. The effect on the total (iso)spin j is however 
less clear. Although the action of and of on the states 
give eigenstates of and /a, these are not eigenstates of 

or P. Lengthy calculation using the algebra on 
does show that the acting on a state with j gives two 
states in general: one with j + 1/2 and the other with 
j — 1/2. Nevertheless this can also be deduced from a 
much simpler argument. 

From the algebra in Appendix ^ each Of, carries j = 
1/2 so decomposition of a direct product of a state with 
j with 1/2 gives the direct sum 



These relations together with Eq. (CI) allow us to obtain 
up to an overall phase factor 



^ : /(j±'7^+l)(i±7^ + l),^^ , 
Ctb,m,n,+) ^ ±J '^ 20- + l)(2Vl) ^ ^'^^ 



CA3,m,n^-) - (C6b) 



J ® 1/2 = (j + 1/2) ® (j - 1/2) 



(C2) 



Alternatively a state of total (iso)spin j has a wavefunc- 
tion which is a polynomial of degree / = 2j in (ob + iac). 
Acting with on the state wavefunction will give a 
polynomial of Z + 1 = 2(j + 1/2). This seems to indi- 
cate that a state with j will become one with j + 1/2 
but this is not the com plet e story. Remember that there 
is the constraint Eq. (|cTl), it is therefore always possi- 
ble to combine with one in the wavefunction to 
reduce some terms in the polynomial from degree / to 

1 = 2 (j — 1/2). As a result one must generally have 
both possibilities of j -I- 1/2 and j — 1/2. The action of 

on the states can thus be written as 



at\i,m,n) 



a-t \j,m,n) 



+CAj,m,n,~)\j - i,m± i,n± 



= C'4'''(j,"i,7i,+)|j + i,m± i,n=F I) 



(C3a) 



(C3b) 

where and are numerical functions of (j^m^n). 

To solve for the coefficien ts, o ne us es th e algebra in 
Eq. (B8). For example Eqs. (B8a) and ( B8c ) acting on a 
state give relation connecting Cf and to them selve s 
within a (iso)spin multiplet and Eqs. (BSb) and (B8d) 
relate Cf^ to Cf. Then from the complex conjugation of 
matrix elements 

(j + 95>"^± hn± ^\df\j,m,n)* 

= {j,m,n\af\j + q^,m± ^,n± ^) (C4a) 

{j + q^,m± h\af\j,m,n)* 

^ {j,m,n\af\j + q^,m± ^,nT i) (C4b) 



with q = ± and using Eq. (|C3|), we arrive at the relations 
between the coefficient functions 

Cf{j,m,n,qy = Cfij+q^,m±^,n±i) (C5a) 
Ct{j,m,n,qy = Cj^O' + gi, m ± i, n T ^) (C5b) 



APPENDIX D: ALTERNATIVE FORM OF THE 
SU{2) ANNIHILATION OPERATORS 

SU{2) operators are not used in the final results for 
the coherent states. They are to be replaced by those 
from 50(3). We will show, nevertheless, in this section 
and in the next that there is an alternate form of the 
annihilation operators and that which chosen form of the 
secondary constraint among the various available choices 
discussed in jss) is used does not affect the uniqueness 
of the annihilation operators. This implies in turn that 
there is no ambiguity as to the uniqueness of the coherent 
states that follow from the operators. 

The alternate form has a much closer resemblance to 
those of the simple harmonic oscillator than the unfamil- 
iar form of position operators sandwiched between two 
exponentials. The complexifier is 



C = 



1 

2Xuj 



J' = 



1 

2Acj' 



(Dl) 



The u! is a. fundamental energy scale as discussed in the 
text. The spin and isospin operators are related to the 
rotation generators Ltc by 



li = '2^Lq.i — —SijkLjk) ■ 

The annihilation operators are therefore 

; L- p „ p 



E 

n=0 



(D2a) 
(D2b) 

(D3) 
(D4) 



with [A, B]o = A and [A, B]n+i = [[A, B]„, C] represent- 
ing the repeated application of the commutator n + 1 
times as mentioned in the main text. Let us work out 
the first few commutators in the sum. We begin by using 
the basic commutator in Appendix ^ to deduce that the 
angular momentum operators have the commutators 



[Lbc, ad] = i{5bdac - Scd&b) 



(D5) 
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Then using (see the Appendix of 



J' 



1 

—LhcL 



bcJ^bc , 



(D6) 



one gets 



~ -i{Lcbac + CLcLcb) 

= ^i(^Lcb - ^i{d~ l)Scb]ac 
= i{Ma)b . 



(D7) 



The last hne has been written as a matrix equation so 
that Bstl 



This is true because 

[J^Lbc] = 0. 



(D8) 



(D9) 



The annihilation operators can therefore be expressed in 
terms of A4 



Ab = exp{2i-A^}a 



(DIO) 



Before proceeding further, we will introduce the pa- 
rameter a which was used in to distinguish between 
the three ways that the second class constraint derived 
from Eq. (2.5) can be implemented quantum mechani- 
cally. In brief one can choose among any of the three 
dbTTb = 0, TTbdb = or dbTTb + TTbCib — forms. This 
results in the relation 

TTc — dbLbc + ioidc = Lbcdb + i{a — d + l)dc (Dll) 

between tt;, and Lbc- The corresponding value of a is 
given by a = 0, a = (d—l) and a — (d—l)/2 respectively. 
Using Eq. (Dll), one can work out 



(D12) 



The appearance of -Kb here indicates that there is also the 
need of [-Kb, J^]. The commutator of itb with L is similar 
to that with db 



[Lbc,TTd] — i{5bdT^c — ScdT^b) 



(D13) 



because the L are generators of four-dimensional rota- 
tions. The commutators are 



[^bj^] 



i{LcbT^c + T^cLcb) 

1 

, ~, ,. -I 
2 \ 2 

i{Mn)b 



iiLcb- 7rj(rf- l)'5cfc)7rc 



(D14) 



The last line can be worked out explicitly using Eq. (Dll) 



{M7r)b = -{nbdcTTc - TTc'^cdb + iTTcdcdb) ~ -iid + l)TTb 



+ 



-{tTcTTc + a- d + l)db 



(D15) 



It is favorable to repla ce tTcTTc by which can be 
achieved with Eq. ( Dll ). Contracting this equation by 
TTc from the left gives 



Then 



{Mtt) 



TT^TTc = 4P _a(Q,_ Ci+ l)i . (D16) 

- -^(4j2 + (l-a)(a-d+l))a,, 

+i,(a_i(d+l))^, . (D17) 



The explicit appearance of a indicates that the annihi- 
lation operators depend on which of the three choices 
are chosen as the quantum version of the secondary con- 
straint. This is the potential non-uniqueness of Ab that 
we allured to at the beginning. This ambiguity will be 
addressed in the next subsection. 

Following Q for the purpose of working out A, we let 



then the corresponding M is 

M = 



-(4j2 + (l-a)(a-d+l)) 

1 



(D18) 



It is advantageous to introduce a shift Ai — Ai — il/A. 

Ab = eMsh} i^M2tM}d)b (D19) 

where 



1 i(a-id) 



(D20) 



This form has the convenient property that on squaring, 
it yields a diagonal matrix 



(D21) 



Therefore this together with Eq. (D12) will allow us to 
work out Ab- 

The series expansion of Ab can be most easily cal- 
culated by exploiting the diagonal form of Eq. (D21). 



Clearly the even power terms are all diagonal in the col- 
umn vectors representation of d and tt. The sum of even 
ri power terms will therefore be a series in j7/2Aw multi- 
plying the n — term which is just db 



1 - 
= cosh 



2!(2Aw)2 
J 



4!(2Aw)4 
4!(2Aw)4 



be 



ab 



2\uj 



db 



(D22) 



This even power series is the hyperbolic cosine. The sum 
of odd power terms can be re-expressed as an even power 
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series in A4 by taking one power of A4 to act on sep- 
arately 



M 






2Xlu ~ 


3!(2Atj)3 


^ 5!(2Aa;)5 


1 






2A^ ^ 


3!(2Acj)3 


^ 5!(2Acc;)5 


1 






2X^ ^ 


3!(2Acj)3 


5!(2Acc;)5 



(D23) 



Since 



{Ma)b = jTTfc - i{a - ^d)ab , (D24) 
therefore finally this is 



1 I , 1 



+esL_(«-_d)sinh(— )a, 



This is the closest form that the annihilation operators 
can be reduced to the familiar form of the SHO. As men- 
tioned in Ref. [0, the J dependent coefBcients are pe- 
culiar to compact spaces. 



2. The Second Choice: a = d — I 

The conjugate momenta are now given by 

TTc = Lbcdb = abLbc + i{d ~ l)dc ■ (E3) 

Obviously iTbdb = is satisfied. The form that Ab 
take is 

Ab = e^cosh(— )a, 



Qb 



. ^ 1 . , f J \ 

2 



V2Ac.J"^ 



(E4) 



3. The Third Choice: a = (d - l)/2 

This choice gives a symmetric expression for the 
conjugate momenta 

= -^{abLbc + Lbc&b) = UbLbc + -ji{d - l)dc ■ (E5) 
The annihilation operators are 

f— 

_,e^_Lsinh(^),r, 



(E6) 



APPENDIX E: ARE A^, DEPENDENT ON THE 
CHOICE OF THE SECONDARY CONSTRAINT? 

The explicit a dependence means that this form of Ab 
is dependent on the exact choice of the secondary con- 
straint. Apparently this would lead to a problem because 
it implies that the coherent states are also dependent on 
the choice of a. In fact this is not the case in spite of 
the apparent evidence to the contrary. It is Tib which is 
dependent on a and not Ab- This can easily be seen if 
one substitutes for -kb in Eq. (D25). 



1. The First Choice: a = 

With this choice, the conjugate momenta are re- 
lated to the angular momenta via 



TTc — CLbL 



be 



(El) 



Now rewrite Hb in each case in terms of Lbc and db the 
Ab are identically given by 

Ab = e^cosh(— )a, 



^ d . .( J 

A J \ 2\uo 

sinh ( ^ 



V2aJ«^^- (E7) 



in all cases so the a dependence cancels out. On 
the angular momenta Lbc are more fundamental than TTft. 
The formers are independent on the choice a but not the 
latter |Q. It follows that the annihilation operators are 
independent of the choice of how to implement the second 
class constraint. It follows that the coherent states are 
independent of a. 



This obviously satisfies the constraint dbi^b — 0. 
Then the annihilation operators are 



APPENDIX F: THE ALGEBRA ON SO(3) 



Ab = es-^'- cosh Jflfc 

+.e^^sinh(^)^. (E2) 



The basic commutators of the spin and isospin opera- 
tors with the components of the 50(3) collective coordi- 



nates Rij are 



[Ji,Rjk\ 



"^^ijin Rrak 



(Fla) 
(Fib) 
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After experimenting with the commutators of j^, Ja, 
and /a with , we find it useful to form the following 

= Rii - R22 + i{Ri2 + R21) (F2a) 
= Rii + R22 - i{Ri2 - R21) (F2b) 
= Rii + R22 + i{Ri2 - R21) (F2c) 
i?"" = i?ii-i?22-i(i?i2+i?2i) (F2d) 

(F2e) 
(F2f) 
(F2g) 



= Ri3±iR 
= R3i±iR: 

— -K33 



•-23 



132 



The superscripts are designed with the following algebras 
in mind 







[h,R^^] 


= R^\ 


(F3a) 


[J3,R-^] 




ikR^-] 




(F3b) 


[•kR^"] 




[h,R"^] 




(F3c) 




= 0, 




= 0, 


(F3d) 




= 0, 




= . 


(F3e) 



It is clear that i?^^, R^ and raise the third compo- 
nent of spin by one unit, R , R ^ and R~'^ lower that 
by one unit. Also R'^^, R ^ and raise the third 
component of isospin by one and R" , R'^ and 
lower that by one. and R^^ leave the third compo- 
nent of isospin and spin respectively alone. Finally 
leaves the third component of both spin and isospin un- 
changed. The first and second superscripts show how the 
value of m and n respectively will be modified. Unlike 
the SU{2) operators where a?, change m and n in steps 
of one-half, the Rij change them by one unit at a time. 
The half unit increment or decrement is the source of 
the problem in SU{2) which effectively couples the bo- 
son states with the fermion states under df,. With the 
Am = Art = ±1, one should now be able to cleanly sep- 
arate the fermionic Hilbert space from the bosonic part. 

The remaining commutation relations in terms of I± 
and J± are 





= , 


[U,R^^] 


= 


(F4a) 




= , 




= 


(F4b) 


[J+,i?-±] 


= 2i?"±, 


[i+,R^-] 


= 2i?±o 


(F4c) 


[J-.R^^] 


= -2i?o± , 




= -2i?±o 


(F4d) 




= , 


[U,R^'] 


= -i?±+ 


(F4e) 




= 


[f-,R^'] 


= 


(F4f) 


[J±,R^°] 


= 0, 


[i±,R"^] 


= 


(F4g) 




= T2R°° , 


[It,R°^] 


= T2i?°° 


(F4h) 


[J±,R"'] 


= TR^' , 






(F4i) 


the above algebras are summarized 


in Fig. a 


It shows 



isospin triplet with j = I. Acting with it on a state with 
total spin j is equivalent to a direct product of j'^l which 
can be decomposed to give the direct sum 



j ® 1 = (j + 1) © J ® (j - 1) 



(F5) 



as mentioned in the text. As a result one should obtain 
in general 

RP''\j,m,n) = CP''{j,m,n,+)\j + l,m + pl,n + ql) 
+CP'^ {j,m,n,0)\j,m + pl,n + ql) 
+CP''{j,m,n,~)\j - l,m + pl,n + ql) . 

(F6a) 

Note that the third component of spin and isospin 
are modified according to the specially designed super- 
scripts p, q as mentioned above. The coefficient functions 
can be solved. Similar to those of SU(2), the vanish- 
ing commutators in Eq. (F4) provide recursion relations 
that link C"^ within a (iso)spin multiplet. The connec- 
tions between C^'' and a different ' come from the 



non- vanishing commutators in Eq. (F4). These relations 
severely restrict the possibilities of the coefficient func- 
tions. To actually solve for them, one has to make use 
of the constraints which allow the coefficient functions to 
be equated to actual numbers and not just to each other, 
and also similar relations relating the complex conjuga- 
tion of one coefficient to another as in Eq. (C5) in the 
SU{2) theory. For example by simply forming matrix el- 
ements from Eq. ( p6| ) using the fact that the Hermitian 
conjugate on an RP'^ operator is to change the sign of its 
superscripts 



thus {R++y 

obtain 



R- 



(RPi)^ R-P'-'i (F7) 
^, = R~^ etc, one can easily 



CP«(j,TO,n,r)* ^ C-P'-'^{j+r,m+pl,n+ql,-r) . (F8) 

The details of how to solve for these coefficients are shown 
in Appendix |h[ 



APPENDIX G: ALTERNATIVE FORM OF THE 
SO{3) ANNIHILATION OPERATORS 

In the main text the complexifier expressed in terms of 
the (iso)spin operators are the same no matter whether 
it is the SU{2) or SO{3) collective coordinates are used. 
From the method described there, the annihilation oper- 
ators are 



A,, 



e 



'RijC^^ 



00 ^ ^ 



n=0 



(Gl) 
(G2) 



We will rewrite this in terms of Rij and 11^ to show that 



that each RP"^ is simultaneously a member of a spin and they can be brought into a form closer to the annihilation 
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operators of the SHO. To work out the series expansion, 
we use the basic commutators in Eq. (Fl). Proceeding 
in a similar fashion to the previous SU{2) theory 

= 2i{MR\j . (G3) 

We have again introduced a matrix M . Because com- 
mutes with any of its individual Ji component, the ap- 
plication of the commutator n times can be written as 



[i?,j,j2]„ = (2*)"(X"i?), 



(G4) 



In terms of M. , the annihilation operators take the form 



4 = {eM±^M}R),j . 
Matrix multiplying Rij once with M gives 



iMR)^ 



-{RilUjnl — Rml^il)Rmj — iRi, 



2n, 



The last equality is obtained by using 

[i?ifc, Iljfc] = [Rki,^kj] — iSij 



(G5) 

ij ■ 
(G6) 

(G7) 



which follows from Eq. (2.16) and the constraint |^ 

^ikRjk + Rik^jk = . (G8) 

The commutators with the conjugate momenta are 
more complicated, we need the commutators 



[I A 



Then 



—2i{Jkekv 
2i{Mn)ij 



(G9) 



(GIO) 



where 



Contracting one free index of each of these with them- 
selves lead to 

RilUklRkn^jn = —UuUji 

~^i. Jra^ikra '^^ik^i'Jn^kjn ^ ^^kj^ 

(G14) 



The first equality is possible becauses of Eq. ( |G8D and 
the constraint of RR^^ = 1. It follows that 



fliiUji = jiJiJj - J'^kj) + \i{Jke-i]k + iSij) 



4 i^JiJj J ^i^ 



(G15) 



Proc eedin g similarly for the second equation in 
Eq. ( |G13 ), we also have 



- = \{hij - I^Sij) + \i{ik^^Jk + iSij) . (G16) 

Then after contractin g the above expressions with Rmj, 
the first term of Eq. ( Gil ) is 



IlilU,nlRri 



2 



(G17) 



To work out the se cond term of Eq. ( pll ), we will look at 
each term in Eq. ( |G16 ) and examine the parts in turns. 
First there is the 



Rim{Ik^' 



mjk 



211, 



(G18) 



which can be calculated in a similar fashion to Eq. ( p(j ) . 
Second one encounters 



Rimlm — \imjkRim{RlAlk ~ Rlk^lj) 



(G19) 



From the primary constraints of unit determinant 
we can write 



^mjkRimRlj — ^ilnRnk 



(G20) 



then 



{MIl)ij = —{RilUml — Rml^il)^mj — i^ij 
= —{Ril^ml—^ilRml)^mj 
= —^il^mlRmj — Ril^ml^mj ■ (Gil) 



The constraints Eq. (G8) is used to obtain the last line. 
Let us now attempt to obtain a relation between the 
square of the conjugate momenta and that of the spin 
or isospin operator. Starting with 



\£ijk{Rjlilkl — RkAjl) 



h = \iijk{RlAlk — Rlk^l]) (G12) 



and using the contracted commutation relations Eq. (G7) 
and the constraint Eq. (G8), we can write 



Rilhjl = ^{Jmf^ijrn + iStj) (G13a) 
RuT^lj = ^{Irn^ijrn + iSij) . (G13b) 



(G21) 



This is a useful equations linking J and / via a contrac- 
tion with the coordinate component Rij . This can be 
converted to the reciprocal equation 



Rm.i Jm. — li 



(G22) 



Using these equations and the fact that = J^, all li 
can be converted into Ji 



RirnJm-^j 



Jilj 



JiJm Rrnj 



(G23) 



Therefore the second term becomes 

Ril^7iil^mj — ~^i^JiJm J kirL)Rraj ~t- 2^^ij : (G24) 

so finally 



(Ain)y = i(j,j„ 



(G25) 
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The simplest way to perform the sum of the series in 
the defi-iiition of Ay is to suppress the indices of R and 
n and let 



R = 



n 



(G26) 



| p5[ . Treating them as column vectors allow us to write 
M using Eqs. (|G|) and (^) as 



2 — i 



(G27) 



This form is not particularly convenient so a shift by a 
constant times the identity will be introduced as before 
M=M-i/2 



M = 



1 . 



(G28) 



Now the square of this shifted matrix is diagonal 

M'^=-(^{j'^l-JJ) + \i^l = -j'^l. (G29) 



We have used the compact notation J'^ to represent the 
spin operators on the left. Because of this property, the 
annihilation operators can be rewritten as 



=exp{2ij}(exp{^^}i? 



(G30) 



and can be calculated by separating the even and odd 
power terms as in the SU{2) theory. The even power 
terms will sum to 



1 - 

1 + 
cosh I 







2!(Aw)2 




J' 




2!(Aw)2 


4!(Atj)2 







(G31) 



We have now restored the 5*0(3) indices which have been 
suppressed in the equations above and used Eq. (G29). 



The sum of odd power terms can be written as an even 
power series in A4 multiplying MR 



Therefore the alternative expression of the annihilation 
operators are 



A,.. 



Rk, 



e^Ai^ cosh ( -iLfc, 

WojJik 

'^Tr:^ sinh f 
2 J \\ujJ 

— sinh f 
J \XujJ. 



Rki 



ik 



(G33) 



The J_used here is a shorthand of the operators in 
Eq. (G29) and is only meaningful in the series expan- 



sion in . Aij is now in a more familiar and yet dif- 
ferent form with obvious differences in the J dependent 
coefficients. 



APPENDIX H: SOLVING FOR THE 
COEFFICIENTS FUNCTIONS OF Rtj ON THE 
STATES 



To solve for the coefficients Eq. (F6), we use the com- 
mutation relations between J±, I± and i?^' and the 
S0{3) constraints. The commutators give relations be- 
tween the coefficient functions and the constraints allow 
them to be solved up to some arbitrary overall phase 
factor. 

For example the vanishing commutators provide recur- 
sion relations of the C'"' amongst a (iso)spin multiplet. 
Let us take 



(HI) 



/iM iM^ iM^ 
V " 3!(Atj)3 ^ 5!(At^)5 

/J_ M'^ M* 
VA^ ~ 3!(Atj)3 ^ 5!(At^)5 

Xuj ^ 3!(Acj)3 b\{\Ljf ^ 



1 sinh ( —\ (2ifl - i 



R) 



ik 



R. 

ij 

{MR)k, 

{2in - \k)kj 

(G32) 



and 



then 



-2R 



(H2) 
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J+^+=^|j, m, n) = ^ + m + rl + 2){j - ni + rl - 1) C+=^(j, m, n, r)|j + rl, m + 2, n ± 1) 



i?+±J+|j,m,n) = V(i + m + l)(j - m)i?+±|j, m + 1, 7i) 

+ rn+ — m) C^^{j,ni + l,n,r)\j + rl,ni + 2,n± 1) 



Equating coefficients give 



^y {j + m + rl + 2) {j — m + rl — 1)C^^ (j, m, n, r) — ■\/ (j + m + l)(j — m)C^^(j, to + 1, r) 
for r = +, 0, — wliich ieads to 



C+±(j,r7i+ l,n,+) 
C+±O-,TO + l,n,0) 
m+l,n,~) 



/ j + TO + 3 



J + TO + 1 



C+±(i,TO,n,+) 



'(j+TO + 2)(j-TO-l) 



+ TO + l)(j — to) 

C+=^(j,TO,71,-) 



C+±(j,m,n,0) 



/ j - TO - 2 



r 



(H3) 

(H4) 

(H5a) 
(H5b) 
(H5c) 



Tliese reiates tlie spin muitiplet of j amongst tlie coef- 
ficients with tlie same p, q. Because of their recursive 
nature, one can guess immediately from these relations 
that 



tators ||38| , p is associated with the spin to and q to the 
isospin n. Therefore the q dependence must be with n. 
The symmetry of (J,p) ^ (/, q) further suggests that 



C+^(j,TO,n, +) oc + m + 2){j + m+l) (H6a) 
C+=^(j,m,n,0) cx ^/{j + m+l){i - m) (H6b) 
C+^(j,TO,n, -) cx V(i - - ?7i - 1) . (H6c) 

Note that even though q = ±, the expressions on the 
right do not have any dependence on the q. This depen- 
dence must be in the remaining parts. From the commu- 



C++{j,m,n,+) cx V(i + n + 2)(j + n + 1) (H7a) 
m,n,0) « v/a + « + l)0'-") (H7b) 
C++{j,m,n,-) cx ■ (H7c) 



Combining both the dependence on to and n, it can be 
deduced that 



C+±0- m,n,+) - ^{j + TO + 2)(j + TO + l)(j ±n + 2)(j ± n + l)/(j,p = +, 9 = ±, +) (H8a) 
C+±O',TO,n,0) = V(j + TO + l)(j - TO)(j ± n + 1)0- T = +, 9 = ±, 0) (H8b) 

C+=^(j,m,n,-) ^ - - m- l)(jT«)(j T?^- l)/(j,P = +,9 = ±, -) ■ (H8c) 

The /'s are the functions of proportionality which can only depend on j, p, q and +, 0, — . 

J_i?^^|j, 771, n) = V (j + TO + rl + l)(j — TO + rl) C^^(j, to, n, r)\j + rl, m, 71 ± 1) 

r=+,0,- 

= J_|j, TO, 77) - 27?"±|j, TO,, 77) 

= ^ + - + l)C+^(j, TO - l,77,r) - 2C°=^(j,TO, 77, r)^ \i + rl,TO, 77± 1) .(H9) 

Equating coefficients for r = +, 0, — , we have 

+ 777 + rl + l)(j - TO + rl)C+=^(j, TO, 77, r) = V(j + - + to - l,77,r) - 2C°^(j, to, 77, r) (.HIO) 
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Now using Eq. (H4) this becomes 

(j{2j + rl + 1) - 2m^ C+='= (j, m, n, r) 
For each each value of r this is 



-2^J{j + m + rl + l)(j -m + rl) C°^{j,m,n,r) 



m C^^(j, m, n, 0) 



^/j + mC^ {j,m,n,-) = 



- + 2 C°± (j, m, n, +) 

v/(j + rn+l)(j-m) (j, m, n, 0) 

v/j -rn - 1 C°=^(j,TO,n,-) 



(Hll) 



(H12a) 
(H12b) 
(H12c) 



Therefore the non- vanishing commutators provide relations between (7^*^ with different pairs of superscript p, q. Using 
the expressions above for 



C"^{j,m,n,+) = - V(J + m + 1)(J - + ± " + 2)(j ± n + l)/a +, ±, +) 
CO±(j, m,n,0) = mV(j±n + l)(jT«)/(j,+,±,0) 
C7*'±(j,m,n,-) = v/(i + H(J - ™){J T Tn- l)/(j, +, ±, -) . 

I 



(H13a) 
(H13b) 
(H13c) 



Other similar expressions for C*"^ can be similarly de- 
duced. 

The /'s functions can be determined with the help of 
the constraints. For example with = (3,3), the 

constraints impose 
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0- nO+ 



^33 - R R 



= i . (H14) 



r=+,0,- 



Sandwiching this between bra and ket of \j, m, n) gives 
^ (|C0+(j,m,n,r)|2 + \C'>\j,m,n,T)f 

r=+,0,- 

= E (|C°-(j,m,n,r)|2 + |C°°(j,m,n,r)|2] =1 

(H15) 



since = i?"" and [R°+,R°-] = 0. In general 

this equation contains the squared modulus of six yet-to- 
be determined functions. The existence of the recursion 
relations means that only one function for each r in each 
multiplet needs to be determined. Because of j > n, the 
states 



|j,m,±(j + l)) , |j-l,m,±(j + l)) 



and 



\j - l,m,±j) 

do not exist, it follows that 

C^'±(j,m,±j,0) = C^^U,m,±j,-) 

= C™(j,m,±j,-) = 0. (H16) 

Then one can choose the top state of the j multiplet 
n = j in the first line of Eq. (H15) or the bottom n = —j 
in the second to get 

|C°±0-,m,±j,-|-)|^ + |C°°(j,m,±j,+)|2 

+ |C0°(j,m,±j-,0)|2 = l. (H17) 



This has three functions but not all of them are indepen- 
dent. There exists the commutators 



(H18) 



which relate C°^{j,m,±j,+) to C°"(j,m,±j +) Tak- 
ing steps similar to those in arriving at Eq. ( HIS ) , one 
can get 



C''>ij,m,n,+) = Ti 



ljTn+1 



C°±(j,m,n,+). 



(H19) 



j±n + 2 

This finally gives us 

(2j + 3)|C"°(j,m,±j,+)|2 + |C™(j,m,±j,0)|2 = 1 . 

(H20) 

Now repeating the same using the transpose of the 
(hi) — (3,3) constraint equation 

RI3 + Rls + RI3 = + (i?™)' = i (H21) 

and 



COO(j,m,n,+) = Tjl±Z^c±0(j,m,n,+) (H22) 



j±m + 2 
m 



±m + l)(j Tm) 



(H23) 



to get 



2^' + ^ -|C™(j,m,±j,+)p 



j ^ m + 1 



[j + l)j Tm 



C'^"(j,m,±j,0)|^ = l.(H24) 



With two equations and two unknowns, the modulus of 
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the coefBcients are 



|CO°(j- m,±j- 0)1 = - 



J + 1 



(H25a) 



..„±,- 0)1 ^ V(Zi^±MI^H25b) 



J + 1 



and 

m,±j,+)| 
\C°^U,m,±j,+)\ 



{j + m+l){j - m + 1) 



J + 1 



2j + 3 



l2{j + TO + l)(j - m + 1) 
(2j+3)(j + l) 



(H26) 
(H27) 



These together with Eqs. (H8) and ( |H13| ) arc sufficient for 
the /'s functions to be determined up to a pha se fa ctor. 
For example from the last expression and Eq. (HIS) 



(j + l)v/(2j+3)(2j + l) 



(H28) 



The complete set of coefficient functions with the arbi- 
trary phase factor chosen to be unity are listed in the 
next section. 



APPENDIX I: COEFFICIENTS OF THE ACTION 
OF SO{3) OPERATORS ON THE STATES 

Let us define two common functions of j 



(j + l)v/(2j + l)(2j+3) 



(11) 



and 



1 



jU + 1) 

which will appear in the coefficients below. 



(12) 



CP«(j,TO,n,+) = pq r{])^{j + pm + 2)(j + pm + l)(j + qn + 2){j + qn + 1) (13) 

CP'^{j,m,n,0) = -s{i)y/{j +pm + l)(j ~ pm){j + qn + l)(j - qn) (14) 

C'"^(j,m,n,-) = pq r(j - - pm){j - pm - - qn)(j - qn - 1) (15) 

for p,q = ±. With either p = or q = 

C°%j,m,n,+) = -q r(jV(.? + ™ + l)(i -m + l)(j + qn + 2){j + qn + 1) (16) 

CP°{j,m,n,+) = -p r{j)^y{j+pm + 2){j +pm + + n + + n ~ I) (17) 

C°«(j,m,n,0) - -m.sij)^{j+qn + l)U-qn) (18) 

CP°U,m,n,0) = s(jV(i + pm + l)(j - pm) (19) 

C°''{j,m,n,-) = q r{j - 1)^ (j + m){j ~ m){j - qn){j - qn - 1) (110) 

C''°{j,m,n,-) = p r{j - l)^/{j - pm){j - pm ~ l){j + n){j - n) (111) 

Then with both p ^ q ^ 

C™(i,TO,n,+) = r(jV(j + m + l)(j - m + l)(j +n + l){j - n + 1) (112) 

C°°{j,m,n,0) = -mns{j) (113) 

C™(j,TO,n,-) = rij - l)y/ij + m){j - m){j + n){j - n) (114) 
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